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ABSTRACT 



Torsional shafts made from fibrous composite materials 
are designed for minimum weight. Multiple loading conditions 
are considered. Design variables are the shaft inner diameter, 
the ply thicknesses, and the ply orientations. Design con- 
straints include limits on ply strain, displacement, frequency, 
and Euler and shell buckling. The design task is solved as a 
numerical optimization problem. Examples are presented to 
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INTRODUCTION 



I . 

Optimization and composite materials are two fields that 
have seen considerable developments in recent years. Advance- 
ments in these technologies have provided the design engineer 
with new dimensions, the depths of which are yet to be explored. 

The major material characteristics that are advantageously 
cultivated in composites are the high strength and stiffness 
to weight ratios. The weight reduction is so significant that 
the aircraft and aerospace industries have concentrated major 
research efforts in this field. Investigations in improving 
the resistance of composites to environmental factors such as 
temperature, corrosion, and wear are continuing and findings 
affirm the contention that the future of composites is bright. 
However, other industries are slow in utilizing these materials 
in basic machinery elements, mainly due to a major drawback; 
the high cost of manufacturing and processing. It is in this 
vein that the importance of an optimized design becomes apparent. 
Advancement in fabrication and processing of composite materials 
coupled with the capability to optimize designs may well lead 
composites into a more competitive market. 

The application of numerical optimization techniques to 
the design of composite shafts is demonstrated here. The effects 
of a small mass imbalance on the deflection of a composite 
shaft with synchronous whirl are included. 
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Section 2 defines the scope and limitations of the analysis 
and optimization. 

Section 3 describes the analysis of composite shafts, 
including : 

A. Composite properties 

B. Shaft section properties 

C. Composite failure criteria, ply strain calculations 

D. Displacement and vibration analysis 

E. Summary of failure modes considered 

Section 4 presents the optimization results [Refs. 1, 2]. 

Three design examples are shown. 

Section 5 contains the conclusions drawn from this study 
and recommendations for future investigations. 

Appendix A contains a description of the analysis portion 
of the FORTRAN program used for composite drive shaft design. 
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II. SCOPE AND LIMITATIONS 



The analysis presented here applies to the shafts made of 
isotropic or composite materials with p inned-p inned end con- 
ditions. Basic shaft design formulas which include the effects 
of a small mass imbalance are used. The shaft is assumed to 
rotate with synchronous whirl. 

The optimization is constrained to subcritical speeds only. 
Supercritical speeds introduce feasible design regions that 
are disjoint from the primary design space. In this investi- 
gation, these secondary areas are considered infeasible. 

Since the optimization tends to a large radius-small thick- 
ness design, radial stresses are neglected. However, buckling 
as a thin cylinder due to torsion or compression is considered. 
The equations used for these failure criteria do not include 
dynamic effects. A constant torque, centrifugal force, and 
axial force are considered in the formulas [Ref. 3]. 

The analysis incorporates the ability to design a shaft 
that can be used in two or more loading conditions. For 

example, a shaft may be designed to transmit 50 HP at 300 RPM 
with an axial load of 2000 lbs. as well as to transmit 150 HP 
at 3000 RPM with a 1000 lbs. axial load. 
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III. ANALYSIS 



The problem considered here is a shaft required to transmit 
a specified horsepower at a given speed. Basic equations are 
used throughout the analysis. 

The inch-pound-second system of units (IPS) is used in the 
computer program analysis. However, by changing the equation 
for torque and the acceleration due to gravity to reflect the 
International System of Units, the program may be used with 
consistent SI units. 

The assumptions used in the analysis are: 

1. The plys are made of orthotropic material, but need 

not all be the same material. 

2 . S tructur e : 

a. Plane stress in radial direction 

b. Plys are at different orientations and have 
different thicknesses. 

c. Cross-section is constant over entire length 

3. Loads: (see Figure 1) 

a. Axial force 

b. Constant torque 

c. Centrifugal force 

4. Failure modes: 

Typical ply orientation and composite layup are 

shown in Figure 2. 
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Stress element (Figure 2) 





Figure 1. Applied Loads: Axial load (F), 

Mass imbalance (e), Torque (T). 



axis 
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Ply orientation, 9^ 
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igure 2. Typical ply orientation and composite _ a y u p . 



a. Maximum ply strain limits 

b. Maximum displacement 

c. First fundamental frequency 

d. Shell buckling (axial and torsional) 

e. Euler buckling (contained within the displace- 
ment constraint) 

A. COMPOSITE PROPERTIES 

Stress is a measure of internal force within a body. This 
together with elastic constants are the variables for the 
determination of stiffness and strength of a material. The 
mechanism of deformation and failure are interpreted in terms 
of the state of stress and strain. When discussing stress it 
is usually considered to be the average stress over some physi- 
cal dimension. 

In Composites, three levels of average stress are considered 

a. Micromechanical or local stress is that calculation 
based on distinct, continuous phases of fiber, 
matrix and, in some cases, the interphase and voids. 

b. Ply stress is that calculation based on assumed 
homogeneity within each ply or ply group where the 
fiber and matrix are smeared and no longer recog- 
nized as distinct phases. 

c. Laminate stress resultant N or moment M is based on 
an average of ply stresses across the thickness of 

a laminate. The individual plies are smeared. 
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The coefficients, or material constants, of the stress- 
strain relations can be packaged in set of engineering con- 
stants, compliance components or modulus components. 

a. MODULUS is used to calculate the stress from strain 
This is the basic set needed for the stiffness of 
multidirectional laminates. 

b. COMPLIANCE is used to calculate the strain from 
stress. This is the set needed for the calculation 
of engineering constants. 

c. ENGINEERING CONSTANTS are the carryover from the 
conventional materials. Designers often feel more 
comfortable working with the familiar engineering 
cons tants . 

d. There is a direct relationship between the modulus 
and compliance. One is the inverse of the other. 

The stiffness of unidirectional composites is governed by 
the same stress-strain relation that is valid for isotropic 
materials. However the number of independent consitants are 
four for planar composites as compared to two for isotropic 
materials. The stiffness of isotropic materials can be repre- 
sented by the Young's modulus alone because the Poisson's ratio 
is the same all directions. 

For composite materials, Poisson's ratios are not bounded 
and can have very significant effect. Young's modulus alone 
is not sufficient to describe the stiffness of composite 
materials. It is not sufficient for isotropic materials either 
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The ply stiffness related to the global coordinate system 
are [ Ref . 4 ] : 
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and Q , Q , Q , Q : On-axis modulus components referred to 

x x y y xy ss 

the local coordinate system 
Q^j: Off-axis modulus components; i,J = 1,2,6 

B. SHAFT SECTION PROPERTIES 

Axial stiffness, bending stiffness, and torsional stiffness 
are smeared properties. 
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Axial stiffness (AE) 



a . 



NPLY 

AE - 21 «J U V U 

k=l 



( 3 B . 1 ) 



b. Bending stiffness (El) 



NPLY 

EI ■ 21 CQnVk 

k=l 



(33. 2) 



c. Torsional stiffness (GJ) 

NPLY 



GJ = 



k= 1 



«hi>k J k 



(3B.3) 



where 



«iJ 



k 

NPLY 



the component of modulus ; i,j = 1,2,6 

the ply moment of inertia ; k = 1,NPLY 

the ply polar moment of inertia ; k = 1,NPLY 

the ply cross-sectional area ; k = 1,NPLY 

the number of composite plys. 



C. COMPOSITE FAILURE CRITERIA 



For the determination of strength of any material it is 
usual practice to estimate the stress at the time and location 
when failure occurs. In the case of conventional materials it 
is needed only to determine the maximum tensile, compressive, 
or shear stress or the principal stresses and then some obser- 
vation can be made about the failure mechanism. This process 
is relatively straightforward because isotropic materials have 
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no preferential orientation and usually one strength constant 
will suffice. The isotropic material is essentially a one- 
dimensional or one constant material. The Young’s modulus for 
stiffness will suffice because Poisson’s ratio is the same in 
all directions, and uniaxial tensile strength will also suffice 
because the shear strength is taken to be about 50 to 60 percent 
of the tensile strength. 

For composite materials, however, the on-constant approach 
for stiffness or for strength is no longer adequate. Four 
elastic constants are needed for the stiffness of a planer ele- 
ment. Six constants for the strength of unidirectional composites 
are needed. Unidirectional composites have directionally depen- 
dent strengths. The longitudinal strength can be twenty 
times that of the transverse and shear strengths; so for any 
state of stress, all three stress components must be examined 
before a judgement on the cause of failure can be made. The 
specific stress component that is responsible for the failure 
may not be easily identified. Probably all three components 
are responsible. The effect of combined stresses must be sys- 
tematically determined and can be regarded as a way of life 
for composites . 

For composite materials, there is needed a failure criter- 
ian for the unidirectional plies. The strength of a laminated 
composite will be based on the strength of the individual plies 
within a laminate. Successive ply failures are expected as the 
applied load to a laminate increase. The first ply failure 

is followed by other ply failures until the last ply failure 
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which would be the ultimate failure of the laminate. The ply 
stress and ply strain calculations are intended for s trength 
determination. There are two popular approaches for failure 
criteria of unidirectional composites. They are each based on 
the on-axis stress or strain as the basic variable [Ref. 4]. 

1 . The Maximum Stress Criteria 

*- s - ^ 

Failure occurs when one of the qualities is met. 

2 . The Maximum Strain Criteria 

^ x — ^x 

£y — £y 

£ S — ^ 

Failure occurs when one of the equalities is met. 

The maximum stress and strain criteria are not the same. 
Only when Poisson’s ratio of the unidirectional materials is 
zero, the criteria become identical. Conceptually they are 
similar. Each component of stress or strain has its own cri- 
terion and is not affected by the other components. There is 
assumed to be no interaction. 

In this study the ply strain failure criteria was used. 

D. DISPLACEMENT AND VIBRATION ANALYSIS 

In order to do the displacement and vibration analysis the 
deflected shape of beam is assumed as: 
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(see Figure 3) 
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where the beam is assumed to be simply supported. The deflec- 
tion is maximum at x = 1/2 [Ref. 3] 
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(3D.1) 



where : 

K x = f Aw 2 

K 2 = / Ag 

f = mass per unit volume 
A = cross sectional area 

= shaft speed is radians per second 
g = acceleration due to gravity 

e = eccentricity of mass with respect to the axis of 
rotation 

^ = shaf t length 
El = smeared bending stiffness 
F = axial load 

Certain characteristics of this equation are of particular 
interest. If there is no rotation or axial load, the formula 
reduces to the classical equation for maximum beam deflection 
under its own weight. 
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/ _ 5k 2j2 

max 384EI 



(3D. 2) 



If there is no axial load, the equation becomes unstable 

j) 4 

(deflection goes infinity) when El = K^(^) or, solving for 
the critical speed, 



0 = 

C 




(3D. 3) 



which defines the first fundamental frequency of a simply 
supported shaft. 

If there is no rotation, instabilty arises when 



F - ^ E1 
F c “ ~J2~ 



( 3D . 4 ) 



recognizable as Euler's column buckling criteria. 

With both rotation and axial load, the critical speed 
becomes 



c 





EI - F(^) 



(3D. 5) 



To avoid the instability regions, the denominator of 
equation (3D.1) must be given a lower bound greater than zero 
(point A in Figure 4). Additionally, a maximum deflection has 
to be imposed to prevent computer overflow (point 3 in Figure 4). 



These limits have no relevance with the deflection constraint 
used in the design optimization process (point C in Figure 4), 
but are used only to prevent numerical ill-conditioning on the 
computer . 

E. SUMMARY OF FAILURE MODES CONSIDERED 

Failure criteria serve important functions in the design 
and sizing of composite laminates. They should provide a con- 
venient framework or model for mathematical operations. The 
framework should remain the same for different definitions of 
failures, such as the ultimate strength, the proportional limit, 
yielding, endurance limit, or a working stress based on design 
or reliability considerations. Failures in composite materials 
involve many modes, such as fiber failures, matrix failures, 
interfacial failures, delamination, and buckling. Furthermore, 
the various modes interact and can occur concurrently or 
sequentially . 

NLC (6 (NPLY)+4) modes of failure are considered in this 
investigation. The 6(NPLY)(NLC) give the maximum strain limits 
in longitudinal compressive strain, longitudinal tensile strain, 
transverse compressive strain, transverse tensile strain, nega- 
tive shear strain, and positive shear strain. Failure is assumed 
if a ply strain multiplied by a safety factor is greater than 
the corresponding strain limit. 

A limit on the def lec t ion is imposed. A common practice is 
to specify a maximum deflection in inches per foot of shaft 
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Figure 3. Assumed deflected shape of the shaft. 



Deflection 




Figure 4. Limits on the deflection. The solid curve 
represents the deflection values used in 
the analysis. 
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length. For example, as used in this report, 



X = 0.005 inch/foot length 
^ max 



Since the optimization may lead to a thin cylinder design, 
buckling due to torsion and compression are considered. The 
critical torque and critical compressive stress for buckling of 
thin cylinders are [Refs. 5, 6]: 



T 

c 



7M/2 E 
3(1- V 2 ) 0 ’ 75 




^ «hi 

3 





Experimental data indicate that actual failure of the cylinders 
usually occurs below fifty percent of the values calculated 
from these equations. Therefore a large factor of safety is 
needed here. 



Lastly, failure is assumed if a specified percentage of 
the critical speed as calculated by equation (3D. 5) is less than 
the shaft speed. 



The factors of safety used in this study are listed in 



the f o llowing : 

FAILURE MODE 
Strain 

Torsion Buckling 
Shell Buckling 
Frequency 
Deflection 



SAFETY FACTOR 
2 
5 
5 
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IV. OPTIMIZATION RESULTS 



There are a number of optimization programs available, 
each one using different technqiues in locating the desired 
optimum design. COPES/CONMIN [Refs. 1,2] is a versatile pro- 
gram that may be used for sensitivity analysis as well as an 
optimization tool. Provided with a user supplied analysis pro- 
gram (Subroutine ANALIZ) where the objective function, con- 
straints and other relevant parameters are calculated, it 
determines a usable and feasible sector from which a search 
direction is chosen. This choice of search direction is made 
by using such information as the gradients of the objective 
function and gradients of active and violated constraints. This 
iterative process of minimizing/maximizing the objective func- 
tion by changing the design variables is terminated when no 
further improvement can be made. 

In this study, the objective function to be minimized is 
the weight, with the shaft inner diameter and the ply thicknesses 
and orientations as design variables. The analysis subroutine 
incorporates the analysis of composite shafts described in the 
previous sections. Several design examples are shown to demon- 
strate the effects of the different parameters on the optimized 
design. For the examples considered here, the loads are listed 
in Table 1 together with the optimum weight and the constraints 
are listed in Table 2. 
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Table 1: Optimization Examples 



No . 


HP 


RPM 


Axial 

Load 

(lbs) 


Material 


Op t imum 
Weight 
(lbs) 


i 


15(3 


o 

o 

ro 


0 


Steel 


50 . 2 


2 


150 


300 


0 


G/E 


12.5 


3 


1000 


6000 


3000 


Steel 


26 . 5 




1500 


8000 


0 






4 


1000 


6000 


3000 


G/E 


10.8 




1500 


8000 


0 






5 


25000 


110 


40000 


Steel 


24822 . 0 




10000 


30 


5000 






6 


25000 


110 


40000 


G/E 


1872. 5 




10000 


30 


5000 







Examples No. 1 and No. 2 are designs for one load condition. 
Shaft lengths are 120 inches here. 

Examples No. 3 and No. 4 are designs for two load conditions. 
Shaft lengths are 60 Inches here. 

Examples No. 5 and No. 6 are designs for two load conditions. 
Shaft lengths are 240 inches here. 

The results of optimization are presented in Figures 5-16. 



Table 2 : 



Design Constraints 



Numb er 



Cons tr aint 



For NLC = 1 



1 through 6(NPLY) 
6 ( NPLY) +1 
6 (NPLY)+2 
6 (NPLY ) +3 
6 (NPLY) +4 



Maximum Strain Limits 
Buckling due to Torque 
Buckling due to Compression 
Speed 

Deflection 



For NLC - 2 

6 (NPLY ) +5 through 12(NPLY)+4 
12 (NPLY ) +5 
12 (NPLY) +6 
1 2 (NPLY ) +7 
12 (NPLY ) +8 



Maximum Strain Limits 
Buckling due to Torque 
Buckling due to Compression 
Speed 

Deflection 



Etc . for NLC > 2 



28 



COMPCSI'F C^IVFSI-AFT OUTPUT 



NUMBER OF FLYS 
NUMBER PE MATERIAL TYP C S 
NUMBER -,F ICAO CONDITIONS 
ECCENTRIC I TY 



l 

1 

1 

0 .lOOOOF + OO 



DIMENSIONS 

PLY T HICKNFSS P FF CENT DIAMETER ThETA 
INSIDE C.20000c+0i 

l C • 100 COE + 01 IOC. 00 0 • 400 00 c +0 L 0.0 



STIFFNFSS 
AE = 0.3l07lr+09 
El » 0. * B 8 3 8 C + 09 
GJ = O.27187F+09 



LOADS: 

L.C. T ^ 

1 0 . 315 L 3E>05 0.0 



0 .0 



F 



HP ppM 

Ol l 5000 r + 0 3 3. ‘*0000 = +G? 



LOAD CCNDITICN 1 

PLY Epl S.F 

1 0.0 100 . 

CRITICAL S F r cr ) 

MAXIMUM OEFLECTICN = 

WEIGHT 

VOLUME 



EPT 

00 0.0 
0.16278E+03 
Q.2^099F-0i 
0. 31893F+03 
0. 1L310E+04 



S.F. E D L T S.F. 

100. OU 0.23192E-03 5.6L 



Figure 5 



Initial Design for Example 1 



COMPOSITE CRIVESHAFT OUTPUT 



NUMBER OF PI. YS 
NUMBER nF MT^IAL T YP C 5 
NU M 8 F3 nr LQAO CONDITIONS 
ECCENTE I C I TY 



1 

1 

L 

0 • IQ 000 E+ 00 



DIMENSl ?NS 

PLY T^ICK\ C SS dfRCTNT OIAmETE' THETA 
INSIDE G.5o907E*01 

1 0 • 81 78 6 E-Ol 100. GO 0.53543E+01 0.0 



STIFFNESS 
AE - G.^doSA'^+OS 
El = Q.203 7CE+09 
GJ = 0.1^259E+09 



LOADS: 

L • C • 


T 


M 


F 


H.P . 


l 


0 • 3 15 1 3 F+05 0.0 


0.0 


0. 15300£»03 


LOAD COMO IT ION 
PLY 6PL 


L 

S.F. EPT 


S.F. 


=PLT S.F. 


I 0.0 

CRITICAL 


SPEED 


IDO. 00 0.0 

* 0.2971 7E + Q3 


100.00 


0 . 0 46 90E-03 2 



MAXIMUM DEFLECT I ON = 0.7 04-016-02 
WEIGHT = 0.5 01 9 L c +02 

VOLUME = O.L 779SE+03 



. Optimum Design for Example 1. 



Figure 6 



(V O 



COMPCS! r E C^IVrSHAFT U'JTPUT 



NUMBER n F PLYS 
NUM3EH . c v U c «TAL TYP-S 
NUM6E- Or L^AO CONDITIONS 
ECC FnTP l C I TY 



9 

1 

o.iooooe+oo 



dimensions 

PLY TMICKN C SS 




T 0 1 A M 6T6 ° 


i 

X 

n 

> 




INSIDE 


0 .2Q000 P * 0 1 






i C • 2 8 J 0 C E ► GO 


2 6.00 


0.250GGE+01 


0.0 




2 C.25u0GF+O0 


25 .00 


0.3000J ~+0 l 


0.20000P 


♦ 02 


3 Q.25 JOCE+vjO 


25.00 


0. 3 50 0 J F+0 l 


-0. 20300= 


► 02 


4 C.2500CE*0J 


25.00 


0. -*00006+01 


0*900 JOE 


► 02 


STIFFNESS 
AE = 0.1 2 l 93 c *- 09 
Ei = 0.12926^+09 
GJ = 0.37225 t+o3 











10A0S: 

L.C. T 




M 


P 


H.o. 


j r> '« 


1 Q.31513E+C5 


0.0 0 


.3 


0. 1 50 00 p 


► 0 3 0.3000C6 


10 AC CGNOITICN 


I 










PLY FPL 


S.F. 


6PT S 


. F . 


“°LT S. 


F. 


1 0.0 


100.00 


0.0 


100. 00 


J. 1Q582 C — 02 


I 7. 3 9 


2 0. 9 0 8 1 1 c ~ 03 


2L. 00 


-0.9081 l c - J 3 


9j>.I3 


0.972 736-03 


18.02 


3 -O.WoBF-03 


1 8.00 


0.9761 36- n 


9.89 


0. 1L393F-02 


16. 2 r 


4 0. 5 3 1 59 r - G9 


100.00 


-0.531596- J9 


100.00 


-J. 169316*02 


ic. a 7 


CRITICAL SPEED 


= 0. 


20662F+Q3 




MAXIHif. 0EFL3CTICN = 0. 


1 979 L c - 01 








WEIGHT 


= 0. 


633 3^6+02 








VOLUME 


= 0. 


L 13 1 JE + 09 









Figure 7 



Initial Design for Example 2 



COMPOSIT- CPIV C SH4FT njtPUT 



NU M 3 E c 


n F 


PLYS 


4 


NUMBER 




MA^=^[AL T >P r S = 


l 


NUMBER 


° f 


LOAD CCN'MTIi’NS = 


1 


EC CENT C 


rc 


I TV 


0 .10j00E*00 



DIMENSIONS 
PLY THICKNESS 


H=ar.2*IT OU«ETER 


THETA 


INSIDE 


0.) 28J6 c *Cl 




1 0 . 14-J 3 l E-07 


w.C'u 0 . 1 26 OC E+0 l 

50.00 0 . 1 0^2 1 r-*-0 L 


0.0 


2 C.1^072E>00 


0. 32225 T 02 


3 0.U372P + 00 


50. OC 0. 200 35 R+J l 


- >. 32225^0? 


4 0. 1V> 0 1 E— U7 


C.00 0. 20035 c +0 1 


0.90000E02 


stipends 
AF = 0.21 o42 c -»* 38 
El - 0.764tlP+*07 
GJ = 0.6 0 767 F* 07 







LOADS: 

L.C. T 




M 


F 


H.P. 


3 PM 


l U.31513E+05 0 


.0 0 


.0 


0.15000R+03 3.30000 


LOAC CrNO ITICN 
PLY E p L 


1 

S.F . 


EPT S 


.F. 


-DIT S. 


F. 


1 0.3 


ICO. 00 


0.0 


100*00 


0. 3 32 06 p - 32 


5. c > 


2 0.192C7E-G2 


4.46 


-0. 19207E-32 


9.16 


0. 13364E-02 


10. C 2 


3 -0.23425F-J2 


3. 66 


0.2343 5E-02 


2.01 


0.22 406 F- 02 


8. 2 ■ 


4 0.163LIF-06 


IOJ.OO 


-0. 1631 IE-08 


100.00 


-0.51950F-02 


3.5-1- 



CRITICAL S FF c O = 0.1 1 5 2 ^ F + 03 

MAXIMUM OEFLFCTICN = 0.50C04F-01 

WEIGHT = 0.125 30F+Q2 

VOLUME = 0.223 76E+Q3 



. Optimum Design for Example 2. 



Figure 8 



COM PC SI T E 


CRIVFSH4PT OUTPUT 




NUM 9 F" OF 


PLYS s 


1 


NU y SF3 np 


KiTPD|4L T Y P c S = 


T 


NUMBER p F 


LOAD CCNOIT IONS = 


2 


ECCENTP [CITY 


0.10 309 E + GO 



DIMENSIONS 

PLY THICKNESS OFRCEMT 0IA*FT6P T HE T A 
INS IDE C.l 0 J JOE >02 

1 C. 100 00 E+Ql 100.00 0.120 00 £>02 0.0 



STIFFNESS 
AF = 0.113 S 36+ 1 C 
£ I = 0.173746+it 
GJ = 0.12U26+11 



LOAOS: 

L.C. T M 

1 0.10504F>05 0.0 

2 0.113 176+05 0.0 



F H • n • 

0 • 300 006+04 ol 100006+04 

0.0 0.150006+04 



Q DM 

0.60000 
0. 3000 J 



LOAC CTjND I TICN 1 

PLY . FFL S.F. £PT S.F. c *lT S.F. 

I 0.26 32 3 p - 05 ICC. JO 0.0 100.00 J.51823 c -05 100. 0 3 

CRITICAL SP C 6Q = 0.227426+04 

HA X I K*JM DEFLECTION = 0. 105346-01 



LG A C CONDITION 
PLY EFL 

l 0.0 

CRITICAL SFFFQ 
MAXIMA 0EFL6C 
WF IGhT 
V0LUH6 



2 

S.F. 6PT 

100.00 0.0 
= 0.227436 + 04 
TICN = 0. 2G022F-01 
= 0.58471 C + Q3 
= 0. 207356 + 04 



S.F. =PLT -S.F. 

100.00 0.583016-35 100.00 



Figure 9. Initial Design for Example 3. 



>04 
+ 04 
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CQMPCSITF C c lVESHAFT OUTPUT 



NUMBER OF FLYS 
NUMAF9 r F MATE ^ t M TYD=S 
NUMBER " lC ICAO CONDITIONS 
ECCENTRIC I TV 



1 

1 

C .10000 F + CO 



DIMENSIONS 

PLY THICKNESS PERCENT OIAMFTE* 
INSIDE G.99992E+01 
1 0 .995 2 5 E-01 100.00 0.L0096F+02 



THETA 



0.0 



ST l FFN r S S 
AE = 0.5 1 5926+08 
El = 0. 6 5 C 5 9F+ 09 
GJ = 0.95691E+09 



LOADS: 

L.C. T M 

1 0.1050 9E + 05 0.0 

2 0.118176+05 0.0 



F 

0 • 30Q006 + Q V 

0.0 



H.o. 

0. 100006+09 
0. I 50Q0 C + 04 



RPM 

3.60000 r *19 
C. 3 0000 r + 09 



LOAD CONDITION 
PLY -PL 

1 0 . 5 d 2 C 5 P - 09 

CRITICAL SPEED 
MAXIMUM DEFLECTION 



S.F. EPT 

17.13 0.0 

= 0.2 C6 73 E «-09 
= G. 1 30396-01 



S.F. =PLT S.f. 

100. 00 0 • L 1696E*03 11.16 
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LOAD CONDITION 2 

PLY E PL S.l 

1 0.0 ICO 

CRITICAL speed 
MAXIMUM DEFLECTION = 
WEIGHT- 
VOLUME 



EPT S.F. 

00 0.0 100.00 
0.20691 F + Q9 
0.250072-01 
0.269596+02 
Q.93807E+Q2 



c PL T S.F. 

0.131026—03 9.92 



Figure 10. Optimum Design for Example 3. 
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COM PC SHE CPIVESHAFT OUTPUT 



NUMBF= CP PlYS 
NUMBFR C c M A x 2 F f A L ~Y P-3 
NUMB F R C’P ICAO CCNOITIPMS 
ECCENTP IC I TY 



4 

1 

2 

G.iOOOOE+OO 



DIMENSIONS 
PLY TniCKfirSS 

INS I 0 E 

1 C.2500CF+00 

2 0.250006+00 

3 0.250 jCF +00 

4 C.2500CF+00 



3fRC r M 01 AM ETcR 
C.oO'JOOF+0 l 
25.00 0 • e 5000 -+0 I 

25.00 C.700CC r +Gl 
25.00 0.75000 -+91 

2 5.00 0. oOQQO c + J L 



^HEIA 



0.0 

0.20OC0=+02 
- 0 . 20000 c + 0 2 
0. 90000-^02 



STIFFNESS 
A 6 - 0.298 5 5E+ 09 
FI = O.W21I r +lQ 
GJ = C.4:>C6e c +C9 



LOADS: 

L.C . 

1 

2 



T M 

0.10504F+05 0.0 

0.U817F+05 0.0 



f 

0 .30000 E + O-r 

0 .0 



H P 

ol L00 00 c +04 
0. 150CJP+04 



3 D\| 

). 60000 - + 04 
C. 0000 " +04 



LOAD CCNO I TIGN 
PLY FFL 

1 C. 1 0 C 49 c - 04 100. CC 

2 0.35C6C—C4 IC3.0C 

3 -0.192C5F-04 100.00 

4 0.29259=-I0 UO.Oj 

CRITICAL SP2 C 0 



1 

^.F. »=PT S.c # 

0.0 100.00 
0. 35 03 1S-04 I TO. 00 
0. 2925 4 C — J4 100. JO 
0. 10C49F- J4 100.00 
= 0.201 3UE+04 



MAXIMUM OEFLECTICN = Q.13325E-01 



~°LT S.p. 

3 • 7 5 71 6F- J4 l 0 0. CO 
J • 5 6 0 0 5 c - l>4 100.0 0 
0. 733 84E-0* 100. o; 

-0.9318«E-04 100. C., 



LOAD CCNO f TICN 



2 

S.F. 

100. JO 



PLY E PL 

1 0.0 

2 0.2 9482E- 04 100.00 

3 -0.3 15 6 6 C - 04 ICO. 00 

4 0.32917 c - 10 

CRITICAL SFECO 

MAXIMUM DEFLECTION = 0.2oo95 c -01 
WEIGHT = 0.7 38 90 F +02 

VOLUME = 0.13195E+04 



EPT $ 

0.0 

-0.294628-04 
0 . 31 5c aF- 
1 DU • 00 -0.329L 7E-10 

= J.201 34E+04 



. 8 . 

100.00 

100.00 

100.00 

100.00 



E °t T S.F. 

J. 85181 F— 04 100.00 
J • 702 72 E— 04 100.0J 
0.75291 E-04 100. O'- 
-0.104848-03 iOO.Cw 



Figure 11 



Initial Design for Example 4 



COMPOSITE CRTVESHAFT OUTPUT 



NUMBS 0 r J c °L YS 
MJM0 jF MATE 0 I AL Ty?~s 
NUMBER ■*= LOAD C CNQ I T I ONS 
E.CCENTPICI TY 



4 

1 

2 

0 .lOOOQP+OO 



DIMENSIONS 

PLY THICKN C $S PFRC-NT D! AM PTE° 
INSIDE 0.544318+01 

1 0.15J94E+00 e?.d7 0.5T450E+0’ 

2 0 • 15b 0 6 F-Ol 6.57 0.87762 ~+0 L 

3 0 • 15c J 6 P-Gl d • 5 7 0. 58074^+01 

4 0.592d7F-U 0.00 0.590748+01 



THETA 

0.0 

J . 964Q0P + 0 i 
-0.364008+01 
0. 9QQC0F+ 02 



ST I F FN C S S 
A E = 0.6 7 152P+G8 
El = 0. 2 6 6 5 dE+ 09 
GJ = Q.ld5S9S+G3 



LOADS: 

L.C. 


- 


M 


F 


H.P. 


3PM 


1 


0.1C5CKP+05 


0.0 


0.300 0QE+Q4 


0. 100008+04 


0.60000 r +04 


2 


0.118175+05 


J.O 


0.0 


0. i5000 c +04 


:. 3oooo c +04 



LOAC CONDITION 



PLY 

1 

2 

3 

4 



FPL 

0. 44 542 c - 04 
0 .2 8 5 EOF- 0 3 
— 0. 2 0 0 C 3 C - 03 
0.5 I4S1 C -0S 



C°rTICAl SPEED 



l 

S.F. 

ICO. 00 
29.99 
92. b4 
100.00 



C PT 3 

0.0 

-0.2412 5 P-0 3 
0.2445 75- J i 
0.4454 2E- 04 



- 0 . 2 06 7 1 c + 04 



MAXIMUM 0 E PLEC T t CN = 0.13J3dE-0l 



.F. 

100.00 
72.95 
19. 2b 
100.00 



- PL T S.F. 

3.16223F-0 2 11.34 

0.1 54*t 3 C - 02 11.91 

0 . 15 792E-02 11. 69 

-0. 16400E-02 1 1.22 
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LCAC CONO I T ION 



PLY FFL 

1 0.0 

2 0.2 7 2 54E- 03 

3 -0.2 74C2 c -03 

4 0. 5 79 2 8F- 09 

CRITICAL SFF C Q _ 

MAXIMUM OFFLPCTICN = 0.24942 c -0l 
WE IGHT = 0.108 UE + 02 

VQLUMP = 0.193 14P+03 



9 

S.F. EPT 

ICO. 00 0.0 

31.44 -0.2725 4E- 03 

31.28 0.2740 2E- 03 

100. 00 -0.5792 8E-u9 

- 0.2071 3E+04 



S.F. 

100.00 
6**. 5 d 
17.19 
100.00 



c P L T S.F. 

0.1925 IF-02 10. 03 

Q. 1 7 52 2 --02 1 0. 50 

0 . 17 6 1 7 C - 02 10.44 

-0 . I 9 450E-02 9.97 



Figure 12. Optimum Design for Example 4. 
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COMPOSITE CP l V ES HA F T OUTPUT 



NUMB PR 3P PLYS 
NUMBER 0<= M&TcpiaL TY P r $ 
NUMBER 0^ L n/\Q CONDITIONS 
ECCENTP ICITY 



L 

1 

2 

C .100000 + 00 



DIMENSIONS 

PLY THICKNESS P EP CENT OIAmcjep THETA 
INSIDE C.5000Cc+OL 

1 0 .800 00 E+01 100.00 0.210CCE+02 0.0 



STIFFNESS 
AE = 0.10771 F+ 11 
El - 0.3137 l p + 1 2 
GJ = 0.21 9 tCE+ 12 



LOADS: 

L. C . 

1 

2 



T M 

0<.143?4E + G8 0.0 

0.21008 E+ 08 0.0 



F 

0 .400 C0 C + 05 
0.50000E+0V 



H P 

oT ? 5000F ♦0 5 
0. 10030E+05 



PPM 

0.11000"03 
0. 3 3030 c + 02 



LOAO CCNOITIGN l 

PLY r PI '.F. EPT 

1 0.3Tl36=-05 100.00 0.0 

CRITICAL SFF=0 = 0.1<*6 3<-E*03 

MAX I M*JM ' OE FLECTI CN = 0.1317SF-01 



S.F. r PLT S.F. 

100.00 0 . 6 84905-03 l.R'J 



LOAC COND I T ION 2 
PLY FFL S.F. C PT 

1 0.4642CF-06 1C0.00 0.0 

CRITICAL SP ccr > = 0 . L 96 4 3F +03 

MAXIMUM DEFLECTION = 0.127256-01 
WEIGHT = 0. 2 21 1 3E + 05 

VO LUMP = 0.784146+05 



S.F. EPLT S.F. 

100.00 0 . 10045 6-02 1.2 * 



Figure 13 



Initial Design for Example 5 



COMPOSITE CPTVF SHAFT OUTPUT 



NUMB FP OF PLYS 
NUMB c P Of V ATE - l At. ^Y P C S 
NUMBER 0 p LOAO CCNOITI n N$ 
5CC5NTF I C l TY 



l 

1 

2 

0 . 10 UOO F *00 



DIMENSIONS 

PLY T H I C KNES S PPRC P NT DIAMETF^ THETA 

INSIDE 0. 120005 + 02 

l 0.635395+01 100.00 0.2^7185+02 0.0 



STIFFNESS 
AE = 0 .12091 F+ 11 
FI = 0.57C52E+1 2 
GJ = 0.39 9 265+ 12 



LOADS: 

L . C . 


T 


N 


F 


H • 3 • 


'PM 


1 


0.143 24F+ C3 


0.0 


0.^00 00 = + 0‘5 


0 . ? 50005+05 


0.1 1000 r +03 


2 


C.210095+C8 


0.0 


0 . 500 OOF*Q<* 


0. 1 0000 c +05 


3. 30000 5+02 



LOAD CONDITION 1 
PLY FPL S.F. EPT 

l 0.33083 5-05 1.00.00 0.0 

CRITICAL SPE=0 = 0 .249 93 F +03 

MAXIMUM. 06FL5CTICN = 0.312125-02 



S.F. PPLT S.F. 

100.00 0.44328F-03 2.93 



LOAO CONDITION 2 
PLY f*L S.F. C ®T 

1 0.4 1 353F-06 ICO. 00 0.0 

CRITICAL S P C 5D = 0.2500^^ + 0? 

MAXIMUM 0 5 FL C C TI GN = 0.73531^-0? 

WEIGHT = 0.24-3225 + 05 

VOLUME = 0.880225+05 



S.F. <=PLT S.F. 

ICO. 00 0.65014P-03 2. GO 



Figure 14. 



Optimum Design for Example 5. 
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COMPOSITE CSIVESHAFT OUTPUT 



NUM0FO p F FLYS 
NUM8ER Or M ATP RIAL TY P-S 
NUMBER OF LOO CCNDITIDNS 
6CCENTP I C I TY 



1 

0 .10000 E+ 00 



DIMENSIONS 

PLY THICKNESS P C RC. C NT OMMETFR 
INSIDE 0.50000F+DL 

1 C.200 OCE + Ol 25.00 0. <0000 E+0 l 

2 0.200 OCE + Ol 25.00 J.I3000 c + 0? 

3 0.200 00 E+0 1 25.00 0.17000F+02 

4 C.200 OCF+Ol 25.00 0.21000E+02 



theta 

0.0 

0.2 DOOOE-t-02 
- 0.200 ) 0 C +02 
0.90000 c *-02 



STIFFNESS 
AE = 0.33749E+ 10 
El =■ 0.738 56 F* 11 
GJ = 0.27562=+ll 



LOADS: 

L.C. 

1 

2 



T M 

0 . 143 24E + 0 d 0.0 

C.21006E+08 0.0 



F 

0 .400 00 c + 05 
0.50000=*04 



H P • 

ol 2 50 OOE * 0 5 
0.10000E+05 



RPM 

0.1 LOO > c +^3 
0.3000; : + 02 



LOAD CONDITION 



PLY 

L 

2 

3 

4 



E FI 



L 

S.F. 



0. L 0 3 2 3E- 04 ICO. 00 



0.10948-- 02 
-0.14 L J2 



7.8 3 

5.08 



0.L7L34E-08 ICO. 00 



CRITICAL SPEED 



E*T 

0.0 

-0. 1084 55-02 
0. 142105-02 
0.1032 IE- J4 



* 0. 2 20 7 IE + 03 



MAXIMUM DEFLECT I CN = 0.L0422E-0L 



LOO. 00 
16. 23 
3. 31 
100.00 



FOLT S.F. 

0. 2338 7E-02 7. a 7 

0.2581 IE-02 7. ’ J 

0 . 33 906 F— 02 5.^1 

-0. 54569 P- 02 3.7 7 



LOAD CONDITION 



PLY 

1 

2 

3 

4 



2 

S.F. 



E PL 

0.129C3F-05 ICO. 00 
0.LS935E-02 5. 38 

-0.20812F-02 4. 12 

0.25L29F-08 ICO. 00 



FPT S 

0.0 

-0.1592 2 e -Q 7 
0.2082 56-02 
0. 128 78F-Q5 



= 0.22L02E+03 



CRITICAL S FF r D 
MAXIMUM OeFLECTIGN = 0.100595-01 
WEIGHT = 0 • 4 39 L 2 c +04 

VOLUME = 0.78414E+05 



.F. 

100.00 
11.05 
2 . ^ 
100 . 00 



FPL T S.F. 

0.34301P-32 5.36 

0. 37 94 6 P-02 4. ro 

0.496405-02 3.7? 

- 0 • 3 003 5 E- 02 2. 30 



Figure 15 



Initial Design for Example 6 



COMPOSITE CRJV6SHAFT OUTPUT 



NUMBER rf PLYS 
NUMBER r. c w ATP'^IAL T YR C S 
NUMB FR QF LOAD CONDITIONS 
ECCENTRIC l TY 



4 

L 

2 

0 • 10000 E + OO 



DIMENSIONS 

PLY THICKNESS PF«C C NT DIAMETER 
INSIDE 0*1061 1 F * 0 2 

1 0. 1235 iE+OO 3.85 0.10858 -^02 

2 0.154276+01 Ad. 08 0.13Q436+02 

3 0.154? 76+01 4«.J8 0. 17 029 F + 0 2 

A 0.2384 2 F-J6 0.00 0. i 7029 6+02 



THETA 



0.0 

0.422A t c + 02 
-0.422416+02 
0.90Q0QE+02 



STIFFNESS 
AE = 0. I 1 0 loE+ 1 0 
El - 0.2 71C9-+ LI 
GJ * CU3778C6+U 



LOADS: 

L.C. 


T 


M 


F 


H.P. 


RPM 


i 


0.143246+08 


0.0 


0.400 006 + 05 


0.250006 + 05 


0.1 I 003 =*03 


2 


C.210C8E+G8 


0.0 


0 .500006+04 


0. 10000E+05 


0. 30000 E *02 



LO AO CCND I TIEN 



PLY 

1 

2 

3 

4 

CP IT ICAL SPEED 



1 

S.F. 



FPL 

0,363126-04 ICO. 00 
0. 1 3 3 54E- 02 6.42 

-0.15367F-02 5.40 

0.10 1366-08 100.00 



= 0. I 97 5oF + 03 
MAXIMUM DEFLECTION = 0.13017E-01 



FOT S.F. 

0.0 100.00 

-0. 1294 16-02 13.55 

0 . 1623 OF- 02 2.90 

0.3631 IE-04 100.00 



c PL T S.F. 

3 . 2058 36-02 8.94 

0.218055-03 84.33 

0. 3 46 5 9 e — 03 5 3. 09 

-0.322626-02 5.70 



LOAD CONOITICN 



PLY 

1 

2 

3 

4 



FPL 



2 

S.F. 



0.45369 6* 05 IOC. 00 



0. i 9 3 196—02 
-0.23539F- 02 



A . 4 V 
.64 



0. 148666—08 100.00 



EPT 

0.0 

-0. 1927 46—02 
0.23584F-0? 
0.453 7 56-05 



* 0. 198336+03 



CRl TICAL spe c d . . . 

MAXIMUM DEFLATION - 0.124826-01 
WE IGHT - 0. 187256+04 

VOLUME = 0. 3 343 BE'*' 05 



S .F. 

100.00 

9.13 

2.00 

100.00 



EPLT S.F. 

0.301396-02 6.09 

0.36631 c - 03 49.46 

0.45984P-03 40.01 

-0. 473476-02 3.89 



Figure 16. Optimum Design for Example 6. 
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V. 



CONCLUSIONS AND RECOMMENDATIONS 



With all other parameters held constant, an increase in 
speed drives the design to a large inside diameter while 
reducing the thickness considerably. This trend is so strong 
that a reasonable lower bound on the thickness must be given 
to prevent numerical ill-conditioning. 

The introduction of eccentricity lets the optimization 
produce larger values for the design variables. 

The effect of axial load is reflected on the deflection. 

If very small deflections are allowed, the axial load effects 
are reduced . 

The versatility of the program should not be lost in the 
emphasis on composite shaft applications. The objective func- 
tion need not be confined to the weight. If the outer dimen- 
sion is critical to the design, one may set a specific value on 
the diameter and let all other parameters vary during optimization. 

In summary, numerical optimization provides an efficient 
and effective way of creating composite shaft designs. 

The study has shown the feasibility of using numerical 

optimization techniques in the design of composite shafts within 

the limitations imposed in the analysis. Further studies on 

the same design field may be pursued by eliminating some of these 

limitations. For example, the unfeasible region defined by 

the speed constraint extends only to some discrete distance 

and a supercritical feasible design area does exist. It will 
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be interesting to develop optimization techniques to accomodate 
this disjoint design field and to investigate the peculiari- 
ties this will introduce to the anlaysis. 

Another investigation may be concerned with shafts with 
variable thickness and diameter along its length. The inclu- 
sion of radial stresses in the calculations, the use of finite 
element methods in three dimensional stress analysis of the 
problem, and vibration considerations are but some of the aspects 
of research that may be useful in future anlaysis of composite 



materials . 



APPENDIX A 



COMPUTER CODE USAGE 

The FORTRAN program used in this investigation is described 
here. Subroutine ANALIZ provides the basic analysis used in 
the optimization. It reads the initial design descriptions and 
calculates the values of the objective function, constraints, 
and all other parameters necessary to solve the analysis 
problem. COPES/CONMIN updates the design to minimize/maximize 
the objective function, iterating until no further improvement 
in the objective function is possible without violating one of 
the constraints. This process is shown in the flow diagram. 
Figure A-l. 

The global catalog given in Table 3 lists the locations, 
FORTRAN name, mathematical symbol and description of the para- 
meters used in the optimization. These parameters are contained 
in the labeled COMMON block, GLOBCM, which is accessed by 
COPES for optimization. 

There are seven data groups to describe the initial design. 

The first data group contains the title. One card is used here. 
The second data group contains the number of plies, NPLY, the 
number of different materials, NMT , and the number of load con- 
ditions, NLC . One card is used, here. The third data group 
contains the initial inside diameter, shaft length, magnitude 
of mass imbalance, fraction of shaft radius at which the mass 
imbalance acts, and the acceleration of gravity. One card is 



used here. The fourth data group describes the thicknesses 
and the orientations of plies and material type used in the 
analysis. NPLY cards are used here. The fifth data group 
describes the ply longitudinal modulus, EL, ply transverse 
modulus, ET , ply shear modulus, GLT, ply major Poisson's ratio, 
PRLT, and the specific weight, RHO . The ply transverse Poisson's 
ratio, PRTL, is calculated internally. NMT cards are used here. 
The sixth datagroup specifies the ply longitudinal compressive 
strain limit, EPLC, ply longitudinal tensile strain limit, 

EPLT, ply transverse compressive strain limit, EPTC, ply trans- 
verse tensile strainlimit, EPTT, and ply maximum shear strain 
limit, GMLT . The compressive strain limits are negative num- 
bers, the tensile and the maximum shear strain limits are 
positive numbers. NMT cards are used here. The seventh data 
group describes the torque, moment, axial load, horsepower, 
and shaft speed in RPM for each, loading condition. NLC cards 
are used here. 

The horsepower, torque, and shaft speed in RPM are inter- 
dependent values. If two of them are known, the third one is 
calculated internally. The following table summarizes the 
required information for Subroutine ANALIZ. 
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Table 3: Input Data 



Data 



Group 


Information 


Format 


1 


Title 


20A4 


2 


NPLY , NMT, NLC 


8110 


3 


DIAMI, AL, EO, EC, GRAV 


8F10 . 2 


4 


TN ( I ) , THND(I), MTYP(I) 


2F10 . 2 


5 


EL, ET, GLT , PRLT , RHO 


8 F 1 0 . 2 


6 


EPLC , EPLT , EPTC , EPTT , GMLT 


8F10.2 


7 


T, M, F, HP, RPH 


8F10 . 2 



4 0 
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SUBROUTINE ANALIZ (ICALC) 

Flow Diagram 

ICALC = 1 

Call Subroutine SHFT06 

Read Initial Dimensions 
Read Material Properties 
Read Design Conditions 
Print all Input 

ICALC = 2 

Calculate Geometric Parameters 
Call Subroutine SHFT07 

Calculate Objective Function 
Call Subroutine SHFT03 

Calculate Shaft Section Properties 
Call Subroutine SHFT09 

Calculate All Constraints 
Call Subroutine SHFT10 

Calculate Critical Frequency 
Call Subroutine SHFT11 

Calculate Maximum Deflection 

ICALC = 3 

Call Subroutine SHFT08 

Print Analysis Results 



OPTIMIZATION FLOW DIAGRAM 




Figure A-l 
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Table 4: Global Catalog of Parameters 



Global 


FORTRAN 


MATH 




Location 


Name 


Symbol 


Definition 


1 


VOL 


V 


Vo lume 


2 


DIAMI 


d 


Inner Diameter 


3 


WGHT 


w 


Weight 


4-10 






Dummy Storage 


11-30 


TN (20) 


t 


Thickness 


31-50 


THND (20) 


9 


Orientation 

(Degree) 


51-551 


G (500) 


G 


Constraints 
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